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Spherically symmetric, static model of the cosmological voids is constructed in the framework of 
the Tolman-Oppenheimer-Volkov equation with the cosmological constant. Extension of the Tooper 
result (dimensionless form of the TOV equation) is provided for non-zero A. Then, the equation 
is simplified in a —> 0, X —> 0, X/a = const regime, suitable for largest structures in A-dominated 
universe. Voids are treated as an underdensity regions in the static Einstein universe. Both overden- 
sity and underdensity (relative to static universe) solutions exist. They are identified with standard 
astrophysical spherical objects and voids, respectively. Model is tested against observed properties 
(the radius - the central density relation) and density profiles of voids. Analytical formulae for ra- 
dial density contrast profile and radii of the voids are derived. Some consequences for cosmological 
n-body simulations are suggested. Hints on the dark matter/dark energy EOS filling the voids are 
provided. 
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I. INTRODUCTION 

In both Newton's and Einstein's theories of gravita- 
tion a large set of static^ astrophysical objects well ap- 
proximated in spherical symmetry exist: moons, planets, 
brown dwarfs, stars, globular clusters, galactic clusters, 
and so on. Spherical shape of these objects is the most 
obvious and the most common effect of the universal at- 
tractive force. However, existence of roughly spherically 
symmetric |16j underdensity regions in the large-scale 
matter distribution in the universe, named voids [24], 
was an unexpected discovery |15l I31j. We will show, in 
the framework of the non-linearly perturbed static Ein- 
stein universe (thereafter s.E.u.), the existence of static, 
spherically symmetric objects devoid of matter relative to 
the surrounding. They belong to the same family of so- 
lutions as normal overdensity objects ranging from plan- 
ets to galactic clusters. Spherical underdensity regions 
("voids") are however directly related to the existence of 
cosmological constant A. Due to very small numerical 
value of A, they exist only for the largest-scale structures 
of universe: the voids. For smaller objects, A provides 
negligible contribution. With A = only overdensity 
static solutions do exist. 

Despite extensive literature on the static, spheri- 
cally symmetric solutions of Einstein equations (see e.g: 
[H [TBI [22] and references therein) this particular branch 
was not discussed in the literature known to author. Pa- 
pers of [13j and [1 deal with similar model, but only the 
overdensity solutions were found and analyzed. 



*URL: |http: //wMw.rlbes ■ if .uj ■ edu.pl/ > Electronic address: 
|odrzywolek@th.if.uj.edu. pl| 

^ More precisely, quasi-static, as all these objects evolve. 



The article goals are: 

1. to show the existence of the static solutions with 
positive (!) pressure gradient 

2. to derive and analyze basic equations for the gen- 
eralized polytropic fluid spheres with non-zero A; 
this generalize result of Tooper |29J 

3. to emphasize existence (and explore properties) of 
the two classes of the solutions; static Einstein uni- 
verse is a boundary between them 

4. to show repulsive "gravitational force" acting on 
the test particles in the vicinity of the " hole" 

5. to use the above solutions in the spherically sym- 
metric static model of the cosmological voids 

6. to test the model against observed density profiles 
and Pc — Dvoid relation for the voids 

Meanwhile, renewed interest in static Einstein universe 
resulted in a large number of recently published papers, 
see e.g. [lJ[31[5l[71ia[Tni[T21[Tnj[2nillSl[2Sl[3n] and 
references in these papers 

Article is organized as follows. In Sect. [ll] we remind 
well known properties of the s.E.u, whose physical prop- 
erties (density, pressure and radius) are used in next sec- 
tion. Sect. |III| to simplify notation. In the Sect. |III| we 
analyze family of solutions to the Tolman-Oppenheimer- 
Volkov (see e.g. [27j) equation with the cosmological 
constant. Some surprising properties of the underden- 
sity solutions are presented in Sect. |IV[ Generalization 
of the Tooper [29] result (valid only for polytropic EOS) 
is presented in Sect. [V| Finally, in Sect. [Vl] attempt is 
made to fit solutions of the TOV-A (Tooper- A) equation 
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to the observed density profiles of the voids. Possible con- 
sequences of the results are outlined in the last section, 
Sect.Eni 



II. STATIC EINSTEIN UNIVERSE 

Traditionally, static universe is discussed in the context 
of the general Friedmann-Lemaitre-Robertson- Walker 
(FLRW) solution, see [B]. Robertson put it into sep- 
arate class " E" . 

For a brief introduction into s.E.u we restrict ourselves 
to the metric representing 3-sphere with radius R{t): 



ds^ = -dt^ + R{tf [dx^ + sin^ X dn^] 



(1) 



Here, dVi^ = siv? 



dd^, and Xi^^^'P ^re spherical 



coordinates on unit 3-sphere. The Einstein equations are: 



i?2 



= SnGp + A 



(2a) 



\ 87rG(p+3p)=A 
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"overdensity" region , 








Static EiH^tein Universe (p^,Pj;) 


- "hole" region 





FIG. 1: Graphical solution of ([6| with EOS Q. Dashed line 
represents eq. ([6|, and solid line corresponds to EOS. Dash- 
dotted line shows relativistic EOS p = p/3. Static 3-sphere 
solution (line crossing) exist for any "normal" EOS. 



3^ +47rG(p + 3p) = A 

IX 



(2b) 



Additionally, we have to specify properties of mat- 
ter filling the universe. We consider Equation Of State 
(thereafter EOS) of the barotropic ideal fluid: 



(3) 



where p - pressure, p - energy density. 

Some of the examples widely used in cosmology are 
dust {p = 0) and relativistic gas {p — p/3). Later in the 
article we will use sub-class of barotropic EOS, polytropic 
gas: 

p^Kp"', 7=1+1. (4) 
n 

If we assume radius of the 3-sphere constant in ([2]): 
R{t) = Re 

then the Einstein's equations reduce to a simple system 
of algebraic equations: 



The gravitational field "source" p + 3p required for static 
solution does not depend on the EOS. For all possible 
static Einstein universes is exactly the same. We will 
show later, that quantity pE + 3pE is very important for 
inhomogeneous static solutions. 

If we specify standard EOS^ and the value of the cos- 
mological constant A, then solution to ([s]) always exist 
and is unique (cf. Fig. [T]) . The well-known Einstein dust 
solution is: 



with the radius: 



Pe 



Re 



A 

AttG 



In the case of relativistic EOS p ~ p/3 we have: 

A 



pe = 



SttG 



(7) 



(8) 



SttG pe — ^ - A 
Re 



SttG pe = — 



Re' 



■A 



(5a) 



(5b) 



where EOS is given by (|4|. 

Let us note, that linear combination of eqs. ^ gives 
(see also (2bl): 



AttG {pE + ipE) = A. 



Re 



(9) 



Therefore we can see, that size {Re) and total mass of 
the Einstein static universe: 



Me ^PeV = 2tt^Re^Pe 



(10) 



(6) 2 With p,p>0 and p < p/3. 



3 



vary only slightly. Numerically: 



and: 



1/VA <Re < 1.22/\/A 



1.57 ,^ 1.44 

= > Me > ^ 

Gy/A GVA 



(11) 



(12) 



If one allows more "exotic" matter (e.g. with negative 
pressure or density) it is possible to produce solutions 
spanning much more wide parameter range, see (llj . 

It is illustrative to compute Re and pe using the value 
of A recently determined from WMAP and SDSS obser- 
vations. According to [58 cosmological constant is: 



Ac^ 
G 



= 6.95 X lO'^'^kg/m 



and density of the baryonic and dark matter is ri„j+rifc 
0.41 r^A- From (|5l we get: 



Radius of the static model (111 is proportional to 



1 

7x 



14.3 Gpc. 



S.E.u. however, once considered very seriously [8], is 
no longer considered a physical reality. 

In the next section we would like to present alternative 
view of the static Einstein universe: one of the static 
spherically symmetric solutions of the TOV-A equation 
with barotropic fluid. 

Our goal is to show the existence and explore prop- 
erties of the wider than s.E.u. class of static solutions, 
namely spherically symmetric, possibly closed, universes 
with EOS As we will show in the next sections, 

3-sphere homogeneous solution is not an exclusive so- 
lution to the given problem. In particular, we discuss 
Tolman-Oppenheimer-Volkov equation with cosmologi- 
cal constant (next section) , generalized Tooper equation 
with non-zero A (Sect, [v]) a nd additional "void equa- 
tion" ( 28 ) in subsection VI A Set of both numerical and 



analytical solutions to these equations is then used to fit 
observed properties of the cosmological voids. 



III. TOV-A EQUATION 

Using standard static spherically symmetric spacetime 
metric: 

ds^ - -e^''^'-) dt^ + e^^^'-^dr^ + dn^ (13) 

we can easily derive Einstein field equations for static 
ideal fluid with barotropic EOS ([3|. As these calculations 



are extensively discussed in literature and textbooks we 
jump immediately to final form of the equations: 



dp 
dr 



-G 



Aiir 



3 p 



(1 



dm 
dr 



2 Gm 



(14a) 



(14b) 



Eq. (14a I is referred to as TOV-A equation. If we put 
A = in (14a I we get standard TOV equation; if addi- 
tionally we neglect pressure as a source of gravitational 
attraction and drop term 1 — 2m/r in the denominator, 
responsible for spacetime deformation in the presence of 
strong gravitational field, we get standard condition of 
hydrostatic equilibrium in Newtonian theory. 

If we solve system (14) then metric (13 1 components 



are easily obtained from formulae: 



„2A 



1 -2m/r- 1/3 r^A 



dp 
p + p 



Kr) 



(15) 



(16) 



where h is the enthalpy. 

As an exercise, we would like to look for constant so- 
lutions of the TOV-A equation. Let us substitute values 
p{r) = Pe and p{r) = pe into (14l. We get immediately 
condition ([6|. The grr component of the metric can be 
put in the form: 



o2A 



with: 



Re' 



= A - ^TiGp 



E- 



We have got 3-sphere with radius Re equal to value de- 
termined for static Einstein universe (5b), as expected. 



Very interesting question arise: what happen if we take 
an initial value'^ for eq. |l4a : 



p{r = 0) = po ^ pe 



The problem is solvable numerically. However, remov- 
able singularity is present at r = 0. To start numerical 
integration of (14) we follow standard procedure, replac- 



ing unknown functions with a series expansion: 
p{r) c:i po + ar +^(3r'^ + . . . 



(17a) 



^ Initial value for second equation (|14b[l is m(r = 0) = 
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FIG. 2: Some selected numerical solutions to the (14 1 with barotropic EOS Q. Values used to obtain figures are: A — 47r, 
G = 1; respective density and radius of the dust s.E.u are Re ~ and pE ~ 1- Initial values are: po = /Ob/IO (dotted) 

and pb/2 (solid) for a "hole" and 2pE for "overdensity region" (dashed). 



m(r) ~ niQ + Ar + ]^B r'^ + r'^ + 



(17b) 



(17c) 



and: 



Substituting series expansion into ( 14 1 immediately 
gives: 



p{r) ~ po + \pr'' + 



p{r) ~ po + 2^''^ + 



m(r) ~ -7rr po + 



(18a) 



(18b) 



(18c) 



f} = b/ 



dp 
dp 



(19b) 



From ( 19 1 we can distinguish three classes of the solu- 
tions: 

(A) If pe + ipE < Pt) + 3po then 6 < and we get so- 
lution similar to normal "star" (cluster, etc.) type 
solutions where density and pressure decreases out- 
wards. 

(B) Pe + 3pE > po + 3po implies 6 > i.e. both pres- 
sure and density grow with radius. This solutions 
describes "hole" in the background matter. 

Finally, if pe + Sps = Po + 3pa we get static Ein- 
stein universe, boundary solution between "clusters" and 
"voids". 

Because class (B) of the solutions remained unnoticed'* 



where: 



b = -nG{po + po){pE + ipE - Po - 3po) (19a) 



In dimensionless form, eq. (|22|, it is much easier to see the 
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(to my knowledge), comment is required. The TOV-A 
equation usually is viewed as a perturbation to normal 
TOV equation. It is applied to the objects like stars. 
Typical initial value of density at r=0 exceeds that of the 
s.E.u. by many of orders of magnitude. Moreover, the 
integration is finished at point where p = p = 0. There- 
fore, no asymtotic density behavior is observed. Addi- 
tionally, asymptotic behavior depends on polytropic, like 
for Lane-Emden equation. For fractional n solution can- 
not be extendended beyond p < 0. If the initial value of 
source density (p+3p) is however comparable to Pe+-^Pe 
(but still in the regime (A) ), then the new behavior is 
observed. Instead of going to p = 0, density approaches 
constant value, exhibiting oscillatory behavior, cf. Fig.|2] 
It happens for relatively small subset of all possible initial 
conditions. 



Surprisingly, starting the integration of (14) with cen- 
tral density in the regime (B), i.e. less than pE + SpE, 
we get qualitatively different solution. We get family of 
static, underdensity objects, again approaching constant 



density (Fig. 



Analysis of the original TOV-A with 



that the hole repels test particles^ from inside. Without 
positive pressure gradient it cannot be static. This also 
explains why non-trivial EOS, e.g. (4]), is required for 
this type of the solution. Dust particles will be repelled 
from the hole. Static Einstein universe is an exception, 
because of its high symmetry. 



V. GENERALIZED TOOPER FLUID SPHERES 
AND DIMENSIONLESS FORM OF THE TOV-A 
EQUATION 

By generalization of the Tooper [29] result, we can put 
TOV-A equations (14 1 in dimensionless form. This is 
possible because power-law equation of state ^ has been 
used. For general barotropic EOS, eq. (14 1 must be used 
instead. 



Let us introduce new dimensionless radial variable: 



general EOS (3 1 is very difficult task. Therefore, we have g^j^^ jj^g-^ functions' 



concentrated on its dimensionless ( 22 1 and simplified (|28|) 
forms. They are valid for polytropic EOS ^ 



We will show in the Section II VI that the existence of 
these objects is a manifestation of the repulsive gravity 
related to the non-zero value of the cosmological con- 
stant. 



p{r)=p,0{xr, 



p{r) ^pc9{x) 



n+l 



(21a) 



(21b) 



IV. PROPERTIES OF THE VOID MODEL IN 
THE FORM OF THE "HOLE" 

"Holes" (Fig. [2]) are very exotic objects. Density in- 
creases outwards and pressure gradient is positive. New- 
tonian intuition tells us, that matter inside given radius 
should attract particle outside. It is however simple to 
show, that test particle will be repelled from the hole^ 
using relativistic approach. 

Because our solution is spherically symmetric we may 
use the effective potential method: 



Veff = 1 - 



2m 




(20) 



where m(r) and ^{r) are computed from (14b I and (16) 
using numerical solution (Fig. [s]) . If no local extremum 
exist, then circular orbits (both stable and unstable) are 
impossible and "force" is repulsive. This is the case of the 
hole (cf. Fig. [3] right panel). Therefore, we have shown. 



ortui 



existence of the two solution classes, cf. Sect. 
^ Or attracted by the universe outside. Unforuinately, for the 
closed universe there is no clear method to define which part of 
the 3-sphere split by 2-sphere into two regions, is "outside" or 
"inside" 2-sphere. We remind, that this solution is not asymp- 
totically flat, and describe spatially finite space 



m(r) = -j^v{x). 



(21c) 



Here, Pc is the density at the center (r = 0) and Pc is the 
central pressure. 

Now, Eqns. ([T4|) take the form: 



(22a) 



x'^ — 2a X V — X x"^ 



(22b) 



where: 

2 _ 4:TtGpc 



A' 



A 

3^' 



PcC 



,2 ' 



7=1-1- 



Initial conditions for system (22) are: 9(0) = 1 and 
v{0) — 0. If cosmological constant (thus parameter A) 
is equal to zero, and we neglect relativity parameter a, 
system (22) reduces to the Lane-Emden equation. How- 



ever, if A ^ 0, there is no simple Newtonian limit. The 
last term of the numerator in eq. (22a I diverge. We will 
use this property in subsection |VI A| to build the simple 
void model. 



Thus, intriguing question arise about "force" between two holes. 
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FIG. 3: Density profile for overdensity object (upper left) and hole (upper right) with effective potential ( |20[ | with very small 
L below. Stable circular orbit exist very close to the center of overdensity region (lower left), while extremum do not exist 
inside a hole for identical L and E. Stable orbit, however, exist very close to the "surface" of the hole. This behavior is easy 
to understand if we note, that source of gravitational attraction/repulsion is over/under-density relative to values for s.E.u., 
not entire matter inside a given radius. 



Behavior of the solution around x — Q can be explored 
using series expansion. Repeating procedure used previ- 
ously for (14 1, we get: 



e{x) ~ 1 + 



1 + a 



A/a 



(23a) 



Using (24 1 we rewrite (23a) in the form 



9{x) ~ 1 + ^ (A - Ac) x^ 



(25) 



v{x) = -x^ + ... 

Let us note that constant solution 9{x) 
Einstein universe) appears only if: 

A = a(a+ 1/3) = A^. 



(23b) 
1 (static 



Eq. 25 shows that 9{x) decreases only if A < Ac; A = in 
particular. 



The transformation from eqns. (14 1 to (22) clarify the 



number of the independent free parameters in ( 14 1 . Den 



sity profile can be obtained from three-parameter special 



(24) function 9n{x;a,X) using (21al 
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VI. FITTING "HOLE" MODEL TO OBSERVED 
PROFILES OF THE COSMOLOGICAL VOIDS 

Radial density profile of the voids can be determined 
from astronomical observations. We have used data from 
Fig. 4 of [15 to compare static model with observations. 

In general, static model based on Tooper-A equation 
(22 1 provide 4-parameter family of solutions. If EOS is 
given completely, i.e. without any of the free parame- 
ters^, then solution is determined by just two parame- 
ters: central density and A. Solutions are enumerated 
by the central density po and value of the cosmological 
constant A. In principle, value of A is known from e.g. 
A-CDM model. However, real voids are not placed in the 
static Einstein universe, but in A-CDM Universe. Some 
effective Ae// ^ Ac dm could be used instead. Therefore 
we consider A i n (p2| ) a free parameter. Radial profile 
of the void in |14j is normalized arbitrarily in radius. 
Therefore additional free parameter A must be used to 
scale solution of the TOV-A equation. Polytropic EOS 
Q has been used in fitting with additional parameter, 
n= 1/(7-1). 

The most general fit: 



TABLE I: Parameters used in Fig. |4] 

pc n logio a logio A A 

0.028 (fixed) 2 (fixed) -66.38 -65.86 1 (fixed) 

0.028 (fixed) 2 (fixed) -1.3 -0.57 0.65 

0.028 (fixed) 8.6 0.66 2.6 0.027 

0.004 ± 0.015 1.26 ± 0.29 -2.5 ± 1.5 -1.0 ± 0.2 1 (fixed) 



for lower-left figure in Fig. |4] a ^ 1 and A <C 1. Little 
can be told on the polytropic index n. For all four cases 
1 < n < 9. The best fit has been obtained with fixed 
central density (lower- left panel in Fig. |4]). Value of A 
is very large in this case, and polytropic index n > 8. 
This suggest that voids are A-dominated objects filled 
with ideal gas with a lot (2n > 16) of the degrees of free- 
dom. This, unfortunately, contradicts results of subsec- 
tion |VIB[ However, used density profile is the averaged 
one. Profiles for individual voids would work better, but 
there is no method to measure them. Too many degrees 



of freedom in model ( 26 1 combined with limited amount 



of data allow us to provide hints rather than constraints 
on EOS. 



p{r) = pcOn{Ar; a. A) 



(26) 



includes five free parameters: central density of the void, 
Pc, constant n coming from polytropic EOS (Eq. [4]), the 
relativity parameter a, the cosmological constant related 
parameter A and one more parameter. A, used to scale 
solution of the TOV-A equation. 

Representative fits are shown in Fig. |4] Overall qual- 
ity of approximation is good. Four cases are presented. 
Because it is unclear why the data points are missing in 
Fig. 4 of [Hj for r < 0.2Ryoid, we consider two scenarios. 
In first, we use fixed central density. The value is deter- 
mined by extrapolation of data from 0.25 R^oid < r < 
0.85 Ryoid to r = 0. Linear fit gives: 



p{r) = pr + q 



where: 



p = 0.009 ±0.004, g = 0.023 ± 0.003 (la). 

This is consistent (p e (-0.0023,0.021), at 95% con- 
fidence level) with constant density (p = 0) for r < 
0.85Ryoid- Assuming constant density we get: 

Pc = 0.028 ±0.001 (Icr), 



A. The "void equation" 

The case with fixed polytropic index, j — 2 and fixed 
central density (Fig.|4] upper-left) is particularly intrigu- 
ing. Solution is not scaled in radius, i.e A = 1. We have 
just two free parameters: a and A. Best fit parame- 
ters (upper-left panel in Fig. |4| are very small, and are 
poorly determined. Actually, from least-square fits we 
get random pairs of extremely small values. However, 
ratio A/a is almost constant. In this situation we can 
simplify eqns. (22 1 a lot. Assuming that a ^ and 
A ^ and A/a — s- const we get: 



V — px 







(27a) 



(27b) 



where we denoted p = A/a. For n — \ (initial conditions 
the same like for (22)) we get analytical solution for the 
density profile of the void: 



p{r) = p, 



sin X — Sp sin x + 3px 



= Ar. 



and this value has been used in fits. 

Four typical cases are presented in Fig. |4j Upper-right 
panel of Fig. |4] will be examined in greater detail in sub- 
section [VlX] All cases are summarized in Table |T] Except 



For n ^ I we still have to use numerics (see however 
31 and Fig. 5|, but the basic equation is much simpler 



eq 



(compared to (22l) 



+ = 3p = g", i?(0) = l,i?'(0) = 0. (28) 

X 



^ For example, if p = p/3 then a = 1/3 and n ■- 



For q = eq. (28) reduces to the Lane-Emden equa- 
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FIG. 4: Best fit solution of tfie Tooper-A (26 1 and radial profile of void provided by [14! . Points show data extracted from 
"SGP" dataset. 



tion. For q^Q behavior of the solutions® for both x ^ 
and X — > oo changes dramatically. Here we restrict to 
positive solutions only: > 0, g > and n > 0. 
For X ^ 1 we have: 



1 



1 
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and depending on sign of g" — 1 we have constant, increas- 
ing or decreasing solution. For x — s- cx) we have obtained 
asymptotic expression: 



a sinujx 



where: 



(29) 



(30) 



The constant a can be determined from numerical so- 



lution. Expression (29) provides surprisingly good (see 



* Scaling symmetry of the Lane-Emden equa tion C,'''d(x/C), where 
s = 2/(1 — n) is no longer valid if g 7^ in l|28[l. 



Fig. O overall fit to the solution in the entire interval 
[OjOojT Value of the constant a in (29) can be roughly 
determined from initial condition d{W) = 1. This gives 



d{x)^q-{q-l) 



sm uox 



(31) 



Expression (|3ip can be used instead of numerical solu- 



tion for (28 1 if large accuracy is not required. Maximum 



relative error is never larger than 15% for t > Q, and is 
below 7% for 1 < n < 5 and < g < 2. 

The most important in our application (a void model) 
is the leading constant in ( 29 1 . Knowledge of the behav- 
ior at infinity allow us (in contrast to eqns. ([14] [22]) ) to 
introduce the density contrast: 

^-P ^ (32) 

P Poc 



Using: 



we have: 



p{r) = 79(Aa;)", poo = Pc 



i+^.f^^y (33) 

p \ <i J 
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FIG. 5: Example solutions to the "void equation" ( |28[ ) 
(solid line) and approximate analytical formula (31 1 (dashed 
line). Dotted horizontal lines show asymptotic value 
lim^_oo '!?n(a;) = q- 



Au! = 0.86±0.09 and q oo. Because pc oc we have 
obtained void model that is empty at the very center. 
Similar result {n — 3/5, Auj = 0.58±0.02 and q 00) has 
been obtained (Fig. |6] dashed line) with fixed polytropic 
index of n = 3/5. If we fix the central density: 

1/g" = Pc/Poc 

then n ^ 1.05 ± 0.2 and Auj = 0.89 ± 0.09, see Fig. |6] 
dotted line. 

This simplest model do not produce impressive results, 
and little information on the matter filling the void can 
be extracted. Unfortunately, data is missing in crucial 
region r < 0.2Ryoid, as well as for r > 2Ryoid- Ex- 
trapolating data we expect po — 0.028 and po 1 for 
r ^ Rvoid- Sudden change of derivative at r ~ Rvoid 
strongly suggest different EOS, more complicated than 
(II. 



B. Radius-central density relation for the voids 




In addition to the void density contrast profile, D^o 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

r/^void 

FIG. 6: Fitting results obtained with aid of the analytical 
formula (33 311 for the density contrast in the void. See the 



text for the details. 



where d is solution to the eq. ( 28 1 or its approximation 

dsTl). 



We have obtained void model much simpler and con- 
venient compared to original function ( 26 ) which must 
be evaluated numerically. Using (33 1 and (31 1 we can fit 



density contrast using three parameters. Two constants: 
A and w always appear as a product Au, and should be 
treated like a single parameter. 

Typical fitting results are shown in Fig. |6] Quality is 
apparently lower compared to the Tooper-A fits (Fig. |4]). 
But curves in Fig. |6] represent functions with well-defined 
and correct asymptotic behavior. Probably, the poly- 
tropic EOS ^ is too simple to reproduce fine structure 
of the density profile. Obtained values of the polytropic 
index n also represent some average value. 

Three curves in Fig. [6] are presented. Solid line shows 
unconstrained fit. Obtained values are: n = 0.9 ± 0.2, 



Pc relation has been published in |2T]. We can use (31 1 
to derive this relation. We define void radius as a point 
where density contrast becomes equal to 1 for the first 
time. Therefore, problem reduces simply to the first zero 
of the sinbjx/ (ljjx), i.e: 

Auj Rvoid = tt. 
Using uj from (30 1 and l/g" — Spc, we get: 



/nA 



fhA 



Spc'- 



(34) 



Here, Spc is the central density contrast plus 1. Using A 
and 7(=1 -I- I/'t-) as a free fit parameters, we get: 

A = 0.23±0.01, 7 = 2.66±0.07 i.e. n = 0.602±0.025 

According to our model, radius of the void can be ap- 
proximated (see Fig. [?]) using: 



id 



Rq 



with: 



i?o = 17.6 ±0.6, 



0.33 ±0.03: 



Polytropic exponent as large as 7 = 8/3 suggest that 
strongly interacting matter (particles?) fill the voids. 
Ideal gas (7 < 5/3) is excluded in this model, because 
void radius would increase with central density. 



VII. CONCLUSIONS 

We have presented analysis of the spherically symmet- 
ric solutions of the TOV-A equation (14 1. We focus on 
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FIG. 7: Measured radii and central densities of the voids 
(Fig. 8 and Tables 2a-c of Ref. Lindner et al. [H] ) fitted to 
our model equation ( |34[ ). Fit is consistent with n = 3/5 i.e. 
7 = 8/3. 



solution branch representing underdensity regions in the 
static Einstein universe. These solutions exist only if ini- 
tial value obey po + 3po < Pe + '^Pe, where index E refer 
to values for static 3-sphere solution. It is shown that 
underdensity void-like static regions (in the presence of 
non-zero A and non-trivial EOS) are possible in general 
relativity. It contradicts common among cosmologist be- 
lief, that voids can exist only because of the expansion. 
In principle, they could exist even in the static universe 
as well. 

Static solution with density increasing outwards is the 
most striking manifestation of the repulsive gravitational 
interaction in the presence of the cosmological constant. 
These solutions remain notoriously unnoticed and ig- 
nored, to my knowledge. For example, authors of the 
paper [T] seem excited with existence of the static gener- 
alization of the Einstein solution. But in their figures 4-6 
pressure (solid lines) always decrease. This is not true 
for all of the solutions, cf. Fig. [2] We also point out 
that for closed static universe there is no such thing as 
" outward" (" inward" ) direction. In general, however, so- 
lutions of ( p2| ) might not be regular for all values of ti, a 
and A. Singularity, if present, break this "two center" 
symmetry, and can be used to define " outward direction" . 

It is shown that hole repels test particles despite p+ip 
is positive inside. Instead, p + 3p — pE — SpE takes the 
role of gravity source, at least in this particular model. 
It can be either positive or negative. Therefore voids in 
static universe grow not only due to attraction of matter 
by structures outside or by expansion. They can increase 
because test particles are expelled from inside. This, in 



turn, accelerate structure formation outside the void via 
"snow plough" mechanism of jl7j . This behavior cannot 
be correctly approximated using standard N-body simu- 
lation. However, we can possibly mimic repulsive forces 
with small admixture of negative mass quasi-particles in 
N-body codes. These "particles" could represent under- 
density "holes" (lack of the matter) similarly to holes 
in semiconductor physics, representing lack of the elec- 
tron. Structure formation in the static Einstein universe 
is of purely academic interest. Similar effects in the real 
expanding universe are not excluded in the non-linear 
regime, however. 

The static model has been applied directly to the ob- 
served profiles of voids. Despite the real universe is ex- 
panding, model works surprisingly well, cf. Fig. [4] Den- 
sity profile can be fitted even with the simple polytropic 
equation of state in the framework of the Tooper-A equa- 
tion (22 1. This might be interpreted twofold: (1) phys- 



ically motivated toy-model which works even if the real 
object do not satisfy all assumptions (2) manifestation 
of the unknown properties of the dark matter or dark 
energy EOS filling the void and making it static. 

To explore the latter possibility, we have derived sim- 
plified version of the eq. ( 22 1 . This " void equation" ( 28 1 



is well-suited for studying largest cosmological structures. 
It comes from a limit a — + 0, A and X/a — const in 
(22 1 that is neither Newtonian nor A-dominated. Basi- 



cally, ( 28 1 is the Lane-Emden modified by the presence 
of the constant on the right-hand-side. The constant 
q ^ break scaling symmetry present in the Lane-Emden 
equation. Unexpectedly, solution for g > is far less 
complicated than that of Lane-Emden equation. Asymp- 
totic behavior (cf. (29l) is the same for all n's. Approx- 
imate analytical form of the solution (31 1 allow us to 
derive Dyoid — Pc relation (34). Comparing with data 
from void catalogue [21 we get 7 ~ 8/3 in 

We have explored hypothesis, that the voids, are not 
completely empty. We have assumed that they are filled 
with unknown form of the "dark" matter. For simplicity, 
and due to poor knowledge of the dark matter/energy 
physics, polytropic EOS (|4| has been used. Static model 
requires also that "voids" are already decoupled from 
Hubble flow, like e.g. galaxy clusters. By extrapola- 
tion of the Hoyle & Vogeley [TJ result for radial density 
profile to r = 0, and from Dyoid — Pc relation [2T| we 
get polytropic index in the range 0.6 n < 1.25. This 
suggest non-ideal gas, possibly strongly interacting. On 
the contrary, if we use central density as a free parameter 
we obtain (i) nearly empty central region and (ii) much 
larger polytropic index ^ ^ n < 9. This do not exclude 
ideal gas, but its nature, e.g. number of the degrees of 
freedom cannot be determined. 
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